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■ A recently proposed strongly correlated electron system associated to the 

• Temperley Lieb algebra is solved by means of the coordinate Bethe ansatz for 

(N 

Q ■ periodic and closed boundary conditions. 

O " 



I. INTRODUCTION 



Integrable highly correlated electron systems have been attracting increasing interest due 
§ . to their potential applications in condensed matter physics. The prototypical examples of 



such systems are the Hubbard and t- J models as well as their supersymmetric generalizations 
Recently many other correlated electron models have been formulated 0^^J^,|,^||,^ . 
Among these an interesting subclass corresponds to models associated to the Temperley- 
Lieb (TL) algebras [TU|. For such models there exists a well established method to construct 



a series of spin Hamiltonians as representations of the TL algebra and of quantum groups, 
the R matrix associated with the XXZ chain being the simplest example [|TT|]. Later on, 
this approach was generalized by Zhang |T2| to construct graded representations of the TL 
algebra using Lie superalgebras and quantum supergroups. Along these lines a new isotropic 



strongly correlated electron model was obtained [|13|, as well as its anisotropic version with 



periodic and closed boundary conditions [|1J] . In addition, it was shown in that this last 
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choice of boundary generates a quantum group invariant model, in contrast to the traditional 
periodic one. 

Models with quantum group invariance and closed boundary conditions were first dis- 



cussed by Martin |T^ from representations of the Hecke algebra. More recently, by means 
of a generalized algebraic Bethe ansatz, Karowski and Zapletal [|r^ presented a class of 
quantum group invariant n-state vertex models with closed boundary conditions. Within 
the framework of the coordinate Bethe ansatz closed spin chains invariant under Uq{si{2)) 
were investigated by Grosse et al. Jl^. Also, an extension of the algebraic approach to the 



case of graded vertex models |jT8| was analysed in where a Uq{spi{2, 1)) invariant susy 
t-J model with closed boundary conditions was presented. 

In this paper we obtain through the coordinate Bethe ansatz approach the solution of 
the anisotropic, or q-deformed, electronic model proposed in ||14[, for periodic and closed 



boundary conditions. Here, the meaning of closed is that an operator, coupling the first 
and last sites, is introduced into the expression of the Hamiltonian such that we obtain a 
quantum algebra invariant closed system (see for more details). In particular, for the 
closed case the Bethe ansatz equations are derived by extending the systematic procedure 
recently developed in E^ to solve the quantum group invariant closed spin 1 chain associated 



to the TL algebra for the case of a graded vertex model. 

The paper is organized as follows. In section 2 we describe the correlated electron system 
associated with the TL algebra. In section 3 we find through the coordinate Bethe ansatz 
the spectra of the model with usual periodic boundary conditions. In section 4 the Bethe 
ansatz solution is presented for closed boundary conditions. A summary of our main results 
is presented in section 5. 

II. THE MODEL 

The starting point for building the model is a 4-dimensional module V of the Lie super- 
algebra Ug{gi{2/1)) utilized to obtain a representation of the TL algebra. Let be 
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an orthonormal basis of V which carries the following parity, 

[|1)] = [|4)] = [|2)] = [|3)] = 1. (1) 
Everywhere we shall use the graded-tensor product law, defined by, 

(a ® 6)(c ® rf) = (-l)WW (ac ® bd) 

and also the rule, 

(|x)®|y))t = (-l)[l-)l[|J/)l(x|®(y|. 

It is then possible to construct the following unnormalized vector of V ®V 
1^) = g-V2 |4) ^ |i) + qi/2 ^ ^ ^-1/2 ^ _ ^1/2 ^ jg^ 
(^1 = g-V2 (4| ^ + ^1/2 ^ _^ ^-1/2 ^ - ^2| ® (3| 

(2) 

Next, to arrive at a hermitian Hamiltonian we consider the operator 

T = 1^) (^1 

A straightforward calculation shows that 

T' = [2{q + q-')]T 
(T®/)(/®T)(r®/) =T®/ (3) 
(/®T)(T®/)(/®T) =/®T , 

such that T provides a representation of the TL algebra. Following the approach of ref. |TI 



to obtain solutions of the Yang-Baxter equation through the TL algebra, a local Hamiltonian 



can be defined by (see ref. Jl^ for more details). 

Hi^i+i = Tj^j+i, (4) 
where on the N-fold tensor product space we denoted 



Ti,i+i = J®('-^) ® T ® J«(^-^-i) . (5) 

In view of the grading the basis vectors of the module V can be identified with the 
eletronic states as follows 

|1) = I + -) = 4cl|0) , |2) ^ I-) = cl |0) , |3) = 1+) = c+|0) , |4) = |0) 

allowing -ffj^j+i to be expressed as 

Hi^i+i = grii,+ni _(1 - - rii+i _) + - ni^+){l - rij _)r2i+i,+?ii+i _ 

+ - _(1 - 71^+1,+) + gn, _(1 - - rij+i _) 

~ — S^_^_i + Cj'^_,_Cj'^_Cj4.i^_Cj+i_+ + c^i^+Cjtj_i _Cj_Ci+ (6) 

+ gc+_^Q+i,+ni _(1 - rij+i + /i.e. - c+_q+i _ni,+ (l - + /i.e. 

+ Q _c,^i - + h.c. - g~^Q,+c^i_+nj+i _(1 - rij _) + h.c. 

where the c-^'* are spin up or down annihilation (creation) operators, the S'^s spin matrices 
and the n'^s occupation numbers of electrons at lattice site i. This model describes electron 
pair hopping, correlated hopping and generalized spin interactions. In the limit g — * 1 it 



reduces to the isotropic Hamiltonian introduced by Links in [|T3|, which was shown to be 
invariant with respect to gl{2) (g) u{l). 
The global Hamiltonian is given by 

N-l 

H=Y. Hk,k+i + h.t. (7) 

k=l 

where h.t. denotes the boundary term. The usual imposition of periodic boundary conditions 
(PBC),i.e., h.t. = Hn,i, has the effect of breaking the f/q((?/(2))®-u(l) symmetry of the model, 
since H^ i 7^ Hi^n, reflecting the non-cocommutativity of the co-product. However, it was 
shown in [|14[ following |jl5|,|16|,|l^|19[ that for a special choice of the boundary term it is 



in fact possible to recover a quantum algebra invariant Hamiltonian which is in addition 
periodic in a certain sense. We shall call this type of boundary by closed one (CBC) and 
denote it by h.t. = Uq (see section 4 for details). 

In the next sections we will find the spectrum of this Hamiltonian for these two types of 
boundaries (PBC and CBC) through a modified version of the coordinate Bethe ansatz. 



III. BETHE ANSATZ SOLUTION FOR PERIODIC BOUNDARY CONDITIONS 



The case with periodic boundary conditions described by the Hamiltonian 
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H — ^ Tk,k+i, 

k=l 



(8) 



can be mapped into a quantum spin chain of N sites each with spin 3/2. To verify this we 
notice that the local Hamiltoniancan be rewritten as 

/ 

T = 1^) (^1 = 



U 04x12 
Ol2x4 Oi2xl2 



where 



U = 



(14| 
(23| 
(32| 
(41 1 



^ |14) |23) |32) |41) ^ 



q -q 
-q q 



1 1 
-1 -1 



1 -1 q-^ q-^ 



q ^ q ^ 



and the following correspondence has to be understood 



|1> = 



M 




f°1 




M 




(o\ 






■ |2> = 


1 




,|3) = 



1 


,|4) = 





















The spin values are given by the eigenvalues of the operator 



3/2 



1/2 



-1/2 



V 



-3/2 



(9) 



(10) 



and the total spin operator commuting with H is Sj^ 
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N 

I 

k=l 



5^ = J®(^-'=). 



(11) 



Following [^] the spectrum of the above Hamiltonian can be classified in sectors which are 
defined by the eigenvalues of the operator number 



r = -N-S' 
2 ^ 



Let us now start to diagonalize H in every sector, i.e., 



(12) 



In the first sectors r = 0, 1,2 the eigenstates do not move under the action of H, i.e.. 



H ^r=o,i,2 = 0. For this reason they are called impurities [22]. In sector r = 3, we encounter 



the situation where the states \a,k) and \—a,k±l), a & = {^,^, |,|}, occur in 
neighboring pairs. They move under the action of H, i.e., the sector r = 3 contains one free 
pseudoparticle. In general, for a sector r we may have p pseudoparticles and A^i and A^^, 

2 2 

impurities of the type ^ and respectively, such that 



r = 3» + A^i + 2A^-i 

2 2 



(13) 



For the first nontrivial sector r = 3, the correspondent eigenspace is spanned by the 
states \k{—a, a) >= || | ■ ■ ■ |— aa | ■ ■ ■ f > , where k = 1,2, N — 1 and a G S"". We 
seek eigenstates of H which are linear combinations of these vectors. It is very convenient 
to consider the linear combination 



\Q{k)) = k{ 



-3 3 
T' 2 

which is an eigenstate of Uk 



-q 



k{ 



1 -1, 

2' 2 ' 



+ q 



^2' 2 



(14) 



Uk \n{k)) = (g + g-i) \Q{k)) = 2{q + q-') \n{k)) . 



(15) 



From now on we will adopt the convention that Uk 
complex spaces at positions k and A; + 1 
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= Uk,k+i operates in a direct product of 



and also a highest weight state, i.e., S'+\l/ = 0. Moreover, the action of Uk±i on \fl{k)) is 
very simple 

Uk±i \n{k)) = \n{k ± 1)) Uk \n{k ± i)) = \n{k)) 

(16) 

Uk\n{m)) = fc^{m±l,m} 

It should be emphasized that the linear combination (0) affords a considerable simpli- 
fication in the diagonalization of H in comparsion with the traditional calculus employing 
the usual spin basis |2T[]. In fact, we believe that this type of ansatz is quite general and 



could be applied to solve a larger class of Hamiltonians derived from representations of the 
TL algebra. 

We will now start to diagonalize H in the sector r = 3. Let us consider the non-trivial 
case of one free pseudoparticle 

^s = j:Mk)\m)- (17) 

k 

Using the eigenvalue equation H = E^^^^, one can derive a complete set of equations for 
the wavef unctions A{k). 

When the bulk of H acts on |^^(A;)) it sees the reference configuration, except in the 
vicinity of k where we use ( [T3| ) and (|TB|) to get the following eigenvalue equation 

(^3 - Q - Q~^)Aik) = A{k - 1) + A{k + 1) 

2<k<N-2 (18) 

Here we will treat periodic boundary conditions . They demand Un,n+i = Un,i, implying 
A{k+N) = A{k). This permits us to complete the set of equations (|18D for A{k) by including 
the equations for k = 1 and k = N—1. Now we parametrize A{k) by plane wave A{k) = AC,^ 
to get the energy of one free pseudoparticle as: 

= 1 (19) 



Here ^ = e^^, 9 being the momenta determined from the periodic boundary to be ^ = 2-711 /N, 
with / an integer. 

Let us consider the state with one pseudoparticle and one impurity of type |, which hes 
in the sector r = 4. We seek eigenstates in the form 

^4(6,6)= E {Ai{h,k2)\n^{h,h)) + A2{h,h)\n2{h,k2))} (20) 

ki<k2 

We try to build these eigenstates out of translationally invariant products of one pseudopar- 
ticle excitation with parameter ^2 and one impurity with parameter ^1: 



^4(6,6) 



(6)) X ^^3(6) + ^3(6) X 



;(6 



(21) 



Using one-pseudoparticle eigenstate solution (|T^) and comparing this with (^) we get 



^l(^),^2(^,| 



ki(-), koi — , 



k,{^),k2{^,^)) + q 



A:i(^),/^2(|^; 



and 



\n2{k,,k2)) 



k,{^,l),k2il: 



+ 



k^ ( — , -), k2(-' 



ki{^,^),k2{^)) + q 



^i(^,^),^; 



(22) 



A,{k,,k2) = A,ei'e2' 



A2{k,,k2) = A2e2'ei'- 



(23) 



Periodic boundary conditions Ai{k2, N+ki) = ^2(^1, ^2) and Ai{N+ki, N+k2) = Ai{ki, ^2), 
i = 1,2 imply that 



N 



(24) 



When H now acts on \E'4, we will get a set of coupled equations for Ai{ki, k2), i = 1, 2. 
We split the equations into far equations, when the pseudoparticle do not meet the impurity 
and near equations, containing terms when they are neighbors. 



8 



Since the impurity is annihilated by if, the action of H on (pO]) in the case far {i.e., 
(/c2 — ki) > 3), can be writen down directly from (p!8|) : 



(e^-Q- Q-^) A(A;i, k^) = k2-l) + A^{ki, k^ + 1) (25) 

and similar equations for A2{ki,k2). Using the parametrization (^3]), these equations will 
give us the energy eigenvalues 

E, = Q + Q-i + ^2 + ^-' (26) 

To find ,^2 we must consider the near equations. First, we compute the action of H on the 
coupled near states k + 1)) and \0.2{k, k + 2)): 

H k + l)) = (Q + Q-^) |l^i(/e, k + l)) + \n^{k, k + 2))- \n2{k, k + 2)) 

(27) 

The last terms in these equations tell us that a pseudoparticle can propagate past the isolated 
impurity, but in doing so causes a shift in its position by two lattice sites. Substituting (^Tf ) 
into the eigenvalue equation, we get 

(^4 - Q - Q-^) Ai{k, k + l) = Ai{k,k + 2)- A2{k, k + 2) (28) 



These equations, which are not automatically satisfied by the ansatz (^31), are equivalent to 
the conditions 

A^{k,k) = -A2{k,k + 2) (29) 



obtained by subtracting Eq. (p8|) from Eq.(|25|) for fci = A; , = ^ + 1- The conditions (^Of) 
require a modification of the amplitude relation (^): 

^ = -er^ = = -1 or = -1 (30) 

In the sectors 3 < r < 6 we also will find states, which consist of one pseudoparticle 
with parameter ^^.-2 interacting with r — 3 impurities, distributing according to (|I3]), with 
parameters ^j, i = I, 2..., r — 3. 
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The energy of these states is parametrized as in (pG]) and ^^-2 satisfies the condition 
(pop with ^ = ,^1 ■ ■ ■C,r-3 Cr-2- It involves only ^^-2 and ^imp = Ci ■ ■ -^r-s, being therefore 
highly degenerate, i.e. 

This is to be expected due to the irrelevance of the relative distances, up to jumps of two 
positions via exchange with a pseudoparticle. Moreover, these results do not depend on the 
impurity type. 

The sector r = 6 contains, in addition to the cases discussed above, states which consist 
of two interacting pseudoparticles. We seek eigenstates in the form 

*6(ei,6)= E Aih,h)\n{h,k,)) (32) 

fci+l<fc2 

Applying H to the state of (0), we obtain a set of equations for the wavefunctions 
A{ki, k2). When the two pseudoparticles are separated, {k2 — ki>3) these are the following 
far equations: 

(Eg -2Q- 2g-i) A{ki, k2) = A{ki - 1, k^) + A{ki + 1, k2) 

(33) 

+ A{ku k2 - I) + A{ki, k2 + I) 

We already know them to be satisfied, if we parametrize A{ki, ^2) by plane waves (^). The 
corresponding energy eigenvalue is 

Ee = 2Q + + ^1 + + 6 + ^2"' (34) 

The real problem arises of course, when pseudoparticles are neighbors, so that they 
interact and we have no guarantee that the total energy is a sum of single pseudoparticle 
energies. 

Acting of H on the state \Q{k, k + 2)) gives the following set of equations for the near 
states 
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H \n{k, k + 2)) = 2 (Q + Q-i) \n{k, k + 2)) + \n{k -l,k + 2)) 

(35) 

+ \n{k,k + 3)) + Uk+i\n{k,k + 2)) 

Before we substitute this result into the eigenvalue equation, we observe that some new 
states are appearing. In order to incorporate these new states in the eigenvalue problem, we 
define 

Uk+i \n{k, k + 2)) = \Q{k, k + l)) + \Q{k + l,k + 2)) (36) 

Here we underline that we are using the same notation for these new states. Applying H to 
them we obtain 

H \n{k, k + i)) = (Q + Q-1) \n(k, k + i)) + \n(k -i,k + 1)) 

(37) 

+ \n{k,k + 2)) 

Now, we extend (^2]), the definition of \E'6 , to 

^6(6,6)= E A{k,,k2)\n{k,,k,)) (38) 

ki<k2 

Substituting (^) and ( pTl) into the eigenvalue equation, we obtain the following set of near 
equations 

(^6 - Q - Q-^) A{k, k + l)=A{k-l,k + l) + A{k, k + 2) (39) 



Using the same plane wave parametrization for these new wavef unctions, the equation (|39D 
gives us the phase shift produced by the interchange of the two interacting pseudoparticles 

Ai2 i + e + (Q + Q-^) = 6 ^ ^ 

We thus arrive to the (BAE) which fix the values of ^1 and ^2 in the energy equation ( p^ 

e = (66)"^ = 1 (41) 
11 



In a generic sector r with / impurities parametrized by ^1^2 ■ ■ - 6 and p pseudoparticles 
with parameters ^;+i^i+2 ■ ■ ■ C,i+p, the energy is 

Er= {Q + Q-' + ^n + C'} (42) 

n=l+l 

with ^„ determined by the Bethe ansatz equations 

a = / + 1,/ + 2, +p, p>2 
eMe2^---ef = (-l)', P=l 

ef-^^ = (-ir, c = i,2,...,/ 

e^ = l, e = ei6---66+i6+2---W (43) 

The energy eigenvalues and the Bethe equations depend on the deformation parameter g, 
through the relation Q + = 2q + 2q~^. 

IV. BETHE ANSATZ SOLUTION FOR CLOSED BOUNDARY CONDITIONS 



The quantum group invariant closed TL Hamiltonians which can be written as ||14|| : 

N-l 

H=J2Uk+Uo (44) 

fc=i 

where Uk is a Temperley-Lieb operator and Uq is a boundary term defined through of a 
operator G which plays the role of the translation operator 

Wo = Gf/^v-iG-i , G = {Q-U^){Q-U2)■■■{Q~UN-l) (45) 

satisfying [H, G] = and additionally invariance with respect to the quantum algebra. The 
operator G shifts the Uk by one unit GUkG^^ = U^+i and maps Uq into f/i, which manifest 



the translational invariance of H . In this sense the Hamiltonian (Q) is periodic. From the 
physical point of view, this type of models exhibit behavior similar to closed chains with 
twisted boundary conditions (see for example pSf for the case of the XXZ chain). 
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The action of the operator G on the states can be easily computed using 



(p!5|) and (|I6|): It is simple on the bulk and at the left boundary 

G\^{k)) = -Q^-^ \n{k + l)) , l<A;<iV-2 (46) 
but has a non-trivial contribution at the right boundary 

G p{N - 1)) = Y: i-Qr' MN - k)) (47) 

k=l 

Similarly, the action of the operator G~^ = {Q^^ — U^-i) ■ ■ • {Q~^ — Ui) is simple on the 
bulk and at the right boundary 

G-^\^{k)) = -Q-^+^\VL{k-l)) , 2<k<N-l (48) 

and non-trivial at the left boundary 

N-l 

G-' \m) = E i-Q)' Mk)) ■ (49) 

k=l 

Now we proceed the diagonalization of H as was made for the periodic case. As ( PD and 
(0) have the same bulk, i.e., differences arise from the boundary terms only, we will keep all 
results relating to the bulk of the periodic case presented in the previous section. 

Let us consider one free pseudoparticle which lies in the sector r = 3 

7V-1 

VI/3 = 5: Aik) \n{k)) . (50) 

k=l 

The action of the operator U = J2k=i Uk on the states \^{k)) is: 

w|fi(i)) = {Q + Q-^)\n{i)) + \n{2)) 

U \^l{k)) ={Q + Q-^) \^l{k)) + \^{k - 1)) + \^l{k + 1)) 
for 2<k<N -2 
U \VL{N - 1)) = (g + g-i) \^1{N - 1)) + \^1{N - 2)) . (51) 

and using (^61) -(^) one can see that the action of Uq = GUn-iG~^ vanishes on the bulk 

Uo \n{k)) = , 2<k<N-2 (52) 
13 



and has the following contributions at the boundaries 

Af-l Af-1 
k=l k=l 

which are connected by 

Uo \n{N - 1)) = (-Q)-^ Wo \^{\)) . (54) 

Before we substitute these results into the eigenvalue equation, we will define two new 
states 

|O(0)) = Wo , = Wo MN - 1)) (55) 

to include the cases k = and k = N into the definition of \l/3, equation (|50|). Finally, the 
action oi H = U + Uq on the states \^{k)) is 

H |fi(o)) = (Q + g-^) |i](o)) + (-g)^ \n{N - 1)) + |i7(i)) 
H \n{N)) = (g + g-i) |i](iv)) + \n{N - 1)) + (-g)-^ |r](i)) (56) 

Substituting these results into the eigenvalue equation = \&3 we get a complete set 
of eigenvalue equations for the wavefunctions 

Es A{k) = {Q + Q-^)A{k) + A{k - 1) + A{k + 1) 

for 1 < < iV - 1 (57) 

provided the following boundary conditions 

{-QfA{k) = A{N + k) (58) 

are satisfied. 

The plane wave parametrization A{k) = A^'' solves these eigenvalue equations and the 
boundary conditions provided that: 
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E3 = Q + Q-' + ^ + r' 

e = {-Qf (59) 

where ^ = e*^ and 6 being the momentum. 

Let us now consider the sector r = 6, where we can find an eigenstate with two in- 
teracting pseudoparticles. We seek the corresponding eigenfunction as products of single 
pseudoparticles eigenfunctions, i.e. 

^6= E A{k,,k2)\n{h,k2)) (60) 

fcl + l<fe2 

To solve the eigenvalue equation H'^q = Eq'^q, we recall (|TBp to get the action of U 
and Uq on the states \Q{ki,k2)). Here we have to consider four cases: (i) when the two 
pseudoparticles are separated in the bulk, the action of U is 

U \n{ki, k2)) = 2{Q + Q-^) |fi(A;i, k2)) + \n{ki - 1, k2)) + \n{ki + 1, k2)) 

+ \n{k^,k2-i)) + \n{k^,k2 + i)) (61) 

i.e., for ki > 2 and ki + 3 < k2 < N — 2; (ii) when the two pseudoparticles are separated 
but one of them or both are at the boundaries 

u k2)) = 2{Q + Q-') k2)) + \n{2, k2)) + k2 - 1)) 

+ |fi(l,A;2 + l)) (62) 

U \n{ki, N-l)) = 2{Q + Q-^) |^](A;l, N - I)) + \n{ki - 1, AT - i)) 

+ \n{kx + l,N -l)) + \n{ki,N -2)) (63) 

U N -l)) = 2{Q + Q-^) N + |fi(2, N-l)) + - 2)) (64) 

where 2 < ki < N — A and 4 < A;2 < — 2; [iii) when the two pseudoparticles are neighbors 
in the bulk 

U \Vt{k, k + 2)) = 2{Q + Q-^) \n{k, k + 2)) + \n{k -l,k + 2)) + \n{k, k + 3)) 

+Uk+i\n{k,k + 2)) (65) 
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ior 2 < k < N — A and (iv) when the two pseudoparticles are neighbors and at the boundaries 
U 3)) = 2(Q + Q-^) |Q(1, 3)) + 4)) + U2 3)) (66) 

U \n{N -3, N-l)) = 2{Q + Q-^) \n{N -3,N-1)) + \n{N - 4, iV - 1)) 

+UN-2MN -3,N -1)) (67) 

Moreover, the action of Uq does not depend on the pseudoparticles are neither separated 
nor neighbors. It is vanishes in the bulk 

Uo \n{ki, k^)) = for ki ^ I and k2^N -I, (68) 

and different of zero at the boundaries: 

Uo |1^(1, k^)) = - eVq)' ^2)) - {-Q)'''-'Uk, \n{k2 + 1)) (69) 

k=l 

N-l 

- E {-Q)'-'m2,k)) (70) 

fe=fe2+2 

Wo \n{ki, N-l))^ (-g)-^+' Uo |Q(1, k2)) (71) 

where 2 < A;i < - 3 and 3 < A;2 < - 2. 

Following the same procedure of one-pseudoparticle case we again define new states in 
order to have consistency between bulk and boundaries terms 

Uo |n(l, k^)) = |O(0, k2)) , Uo \^{ki, N-l))^ |n(A;i, TV)) 
Uo |fi(l, N-l)) = \Q{0, N-l)) + N)) 

Uk+i \n{k, k + 2)) = \n{k, k + i)) + \n{k + i,k + 2)) (72) 

Acting with H on these new states, we get 

H |O(0, k2)) = 2{Q + Q-^) |O(0, k2)) + |O(0, k2 - 1)) + |O(0, k2 + 1)) 

+ \n{i, k2)) + (-g)^-' \n{k2, N - 1)) (73) 
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H \n{ki, N)) = 2{Q + Q-') \n{h, N)) + \n{h - 1, N)) + \n{ki + 1, N)) 

+ \Q{h, N-l)) + (-Q)-^+2 1^^^^ (74) 

H \n{k, k + l) = {Q + Q-^) \n{k, k + l) + \n{k - l, A; + l) + \n{k, k + 2) (75) 

Substituting these results into the eigenvalue equation, we get the following equations for 
wavefunctions corresponding to the separated pseudoparticles. 

{Ee -2Q- 2Q-^)A{ki, fcs) = A{ki - 1, A^a) + A{ki + 1, fca) 

+A{ki,k2-l) + A{ki,k2 + l) (76) 

i.e., for ki > 1 and ki + 3 < k2 < N — 1. The boundary conditions read now 

A{k2, N + k^) = (-Q)^-M(fci, k^). (77) 

The parametrization for the wavefunctions 

A{k^, k^) = Ai^^i'^e'^ + A^ie^e'^ (78) 

solves the equation (|76]) provided that 

Ee = 2{Q + Q-') + 6 + + 6 + ^2"' (79) 
and the boundary conditions (|77D provided that 

C^ = {-Qf-'^ , = (-g)^-4^ ^ = (-g)^(^-) (80) 

^12 ^21 

where ^ = ^1^2 = e**^^^'^^^-', + 62 being the total momenta. 

Now we include the new states (ff2D into the definition of \E'fi in order to extend (|60f) to 



A{k,M)\^{hM)- (81) 

ki<k2 

Here we have used the same notation for separated and neighboring states. Substituting 
(|65|) and ( |75D into the eigenvalue equation, we get 
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(Ee-Q- Q-^)A(k, k + 1) ^A{k-l,k + l)+ A(k, k + 2) (82) 



(83) 



which gives us the phase shift produced by the interchange of the two pseudoparticles 

A12 i + i + {Q + Q-')ii' 

We thus arrive to the Bethe ansatz equations which fix the values of ^1 and 

I i+c+(g + g-%i' 
= i-Qr""-'^ (84) 

Thus in the sector r = 3p, we expect that the p-pseudoparticle phase shift will be a sum 
of two-pseudoparticle phase shifts and the energy is given by 

Esp^j2{Q + Q~' + ^n + C'} (85) 

n=l 

where 



_ (_r)\N-2p+2 TT j 1 + iaih + (Q + Q ^)Ca 



(6e2---ep)^ = (-Q)^^^-'^+'^ (86) 

The corresponding eigenstates are 

*r(ei,6,-y = E A(A;i,/c2,...,M|f](/ci,A;2,...,M) (87) 

l<A;i<...<fcp<Ar-l 

where |0(A;i, /c2, ^p)) = ®f=i |^(^i)) the wavefunctions satisfy the following boundary 
conditions 

A{k^, k2,..., kp, N + k,)^ {-Qf-^P+^Aiki, k2,..., kp) (88) 

It is not all, in a sector r we may have p pseudoparticles and Ni , impurities of the 
type |, respectively. Since if is a sum of projectors on spin zero, these states are also 
annihilated by Uq . Therefore the impurities play here the same role as in the periodic case. 



It means that for a sector r with / impurities with parameters ^i, and p pseudoparticles 
with parameters ^/+p the energy is given by (p6D, and the Bethe equations do not 

depend on impurity type and are given by 



i+p 



i + Ub + {Q + Q-')L' 



ce^e2--<f = i-in-Qr~''^' n - Tr. T ^^T^-i ; m 



b=l+l,bj^a 

with a = / + 1, / + 2, / + p , p > 1, and 



i + Ub + iQ + Q-')^b 



e'i^i+i ■ -^i+pf-'' = i-iYi-Qr^"-''^'^ (90) 

where ^ = ^1^2 ■ ■ ■ ^i^i+i ■ ■ ■ ^i+p- 

Notice in the BAE (^) the presence of a special "g- term" ((— Q)^"^^'*'^) in comparsion 
with the corresponding ones with usual periodic boundary conditions (^31). In fact, this 
feature also appeared in other models ||16| , p!7| , p!9| and seems to be a peculiarity of quantum 
group invariant closed spin chains. 

V. CONCLUSIONS 

We have applied the coordinate Bethe ansatz to find the spectra of the anisotropic 
correlated electron system associated with the TL algebra. This procedure was carried out 
for periodic and closed boundary conditions and the differences between both cases have 
been remarked. 

We believe that the methods here presented could also be applied to solve a larger 
class of Hamiltonians derived from representations of the graded TL algebra, such as the 



orthosympletic models discussed by Zhang in [12|. This is presently under investigation. 



Another interesting extension of this work would be to adapt the methods employed in 
this paper to solve multiparametric versions of these models. p^j25[] 
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